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I. INTRODUCTION 



Usually, in vacuum, the constitutive relations of classical electrodynamics are D = EqE and 
H = B//io. The electric constant eq (permittivity of free space) alone has no direct meaning in 
4-dimensional spacetime; the analogous is true for the magnetic constant fiQ (permeability of free 
space). However, if we combine both constants, the situation changes. As shown by Post [16 ill . 



for example, it is rather the square root of the quotient of both constants, namely Yq := ^/eo/fXo, 
the vacuum admittance of Yq ~ 1/377 ^7 that represents a scalar in 4 dimensional spacetime 
in arbitrary coordinates; the same is true for its reciprocal, the vacuum impedance (resistance) 
:= 1/Yq. Thus, it is possible to extract 4-dimensional information from both 3-dimensional 
constants, provided they are taken together. 

Moreover, c := l/y/Eofio, the vacuum speed of light, has also a 4-dimensional meaning, even 
though c is only a scalar under Poincare (inhomogeneous Lorentz) transformations. This is obvious 
since in noninertial, that is, accelerated frames c is no longer a constant. In this sense, the vacuum 
admittance has a more fundamental significance than the speed of light. The vacuum admittance 



can be measured by a Weber- Kohlrausch type of experiment, e.g., see Raith [|67l1 and Brown lllOll . 



the speed of light by the well-known methods of Foucault or Fizeau. [19211 respectively (even though, 
strictly speaking, the speed of light is put to a certain constant value in SI since 1983, see [52]). 
Sommerfeld's fine structure constant, the dimensionless coupling constant of the electromag- 



netic interaction, can be written as dl 



where Rk is the quantum Hall resistance (von Klitzing constant) associated with the quantum 



Hall effect. It has been shown 112811 that Rk, like the vacuum resistance Zq, is not influenced by 
the gravitational field. This underlines the fundamental importance of Zq as well as that of Rk- 

Let us now turn to media, namely to dielectric and magnetic media that can be described by a 
local and linear constitutive law. Suppose we consider a spatially isotropic medium. In a frame 
where the medium is at rest, we find D = eeqE and H = B/(yU/io), with the permittivity e and 
the permeability fi of the medium; both, e and /x, are dimensionless and depend in general on the 



frequency of the wave studied ("dispersion"). By the analogous arguments as above, ^/esq/ (/i/io) 
is a 4-scalar in arbitrary coordinates and l/^/eJIe^jM) a speed in inertial coordinates. The absolute 
refractive index n of a medium, see ITD, derives from the latter expression as relation of the vacuum 
speed of light to the speed in the medium. Accordingly, n = Jeil. 



Usually media, in particular crystalline media, behave anisotropically. We could generalize the 
laws above by introducing anisotropic permittivity and permeability tensors Eab and f^ab, respec- 
tively, with a, b, ... = 1, 2, 3. But it is better to start right away with a general local and linear 
constitutive law in the context of a 4-dimensional representation of electrodynamics (Secjn]). This 
guarantees automatically relativistic covariance, and the anisotropic laws mentioned will emerge 
as special cases. 

In this context, it turns out, see (fT5l) . that the corresponding 4-dimensional constitutive tensor 
^XuuK. _ _^uXaK _ _^XuKa j^^^ independent components in general, with A, z/, ... = 0, 1, 2, 3. 
In the center of our present paper is only one component of ■^^'^'"^^ namely its totally antisymmetric 
piece. It is the pseudoscalar a := e^x^an (we sum over all repeated indices), which can be 

formed from the constitutive tensor ^^^'^'^ with the help of the totally antisymmetric Levi-Civita 
symbol e^xuaK- We remind ourselves that 'ex^an = +1 or = —1 depending whether XuaK denotes 



an even or an odd permutation of the numbers 0123, respectively; it is zero otherwise, see [18011 . 
Then, 

5 = ^ (^,0123 _ ^0213 ^ ^0312 _ ^0132 _^ _ _ _ _ ) _ 

Observing the antisymmetries of ^'^'^'^'^^ we find 

s = I ix'"'' + x'''' + x'''' 

+^2301 + ^3102 + ^1203) _ (3) 

The components in the first line are related to the magnetoelectric (ME) effect in an external 

magnetic field, MEb, those in the second line to the magnetoelectric effect in an external electric 
field, MEe. That 5 is a pseudoscalar, indeed, can be seen from the transformation properties of 



the quantities involved, see [|26ll . The dimension of a is \l resistance. 



The constitutive tensor y^^^'^^ was seemingly first introduced by Tamm [|82|] . see also Post [16111 . 



and later discussed by O'Dell [|57(] in the context of magnetoelectric media — that is, media in 
which an electric field E induces a magnetic excitation H and a magnetic field B and electric 
excitation D, see Fiebig [17] for a recent review. All four cited authors assumed a further symme- 



try, namely ■y^^"'^'^ = -y/^f^^" ("vanishing of the skewon part", see [^]). This symmetry emerges, 
as soon as one stipulates that the constitutive relation can be derived from a Lagrangian (thereby 
excluding irreversible processes). Then, in particular, i.e., for y^^^'^'-^ = ^^'^'^^'^^ 

a = 1 (;,0123 + ^0231 + ^0312) _ 



Thus we need only three components of the constitutive tensor for the determination of a in the 
case of the vanishing skewon part. Post argued IgI] (not very convincingly, we should say) that 
the pseudoscalar dH) ought to vanish: 5 = 0. This condition was dubbed "Post constraint" by 
Lakhtakia [|40n . It was not assumed by O' Dell so that he was left with 20+1 independent compo- 



nents of the constitutive tensor, see pVll . p.44. 



Later a fierce dispute arose about the Post constraint. The situation was reviewed by Lakhtakia 



4111 . de Lange and Raab [|l3ll . Raab and Sihvola 



Sihvola and Tretyakov [179 



, and in 



27D, see also [|54, 



62 



, Raab and de Lange n66A . Sihvola 



5511 . where more references to the relevant 



literature can be found. The evidence was mounting that there is no reason to assume the validity 
of the Post constraint in general. This point of view will be shown to be correct in this paper. 

According to a theory of Dzyaloshinskii lll5ll . which was based on an analysis of neutron scat- 
tering data, susceptibility measurements, and symmetry (as derived in 1|9|1), a crystal of Chromium 
oxide Cr203 is the substance par excellence for discovering the magnetoelectric effect. In our 
paper, we consider only single domain crystals. In Secjllll we will describe Dzyaloshinskii theory 
of Cr203 and we will determine the 4-dimensional pseudoscalar a of Cr203 in this framework. 
In SecJTVl after a short introduction on dimensions and units, an overview will be given over 
magnetoelectric experiments with Cr203. Corresponding unpublished measurements by one of 
us (J. -PR.) will be presented in some detail. In Sec|Vl we study the experimentally determined 
magnetoelectric moduli of Cr203. Then, we extract from the data, for the first time, the rela- 
tivistic pseudoscalar a for Cr203. It turns out to be temperature dependent and is of the order of 
5 ~ 10~^/Zo. This is a typical magnitude for magnetoelectric moduli (Borovik- Romanov and 
Grimmer ijs^). Thus, it is small but definitely nonvanishing. This proves experimentally that the 
Post constraint is ruled out as a generally valid law. In SeclVTl we show that the pseudoscalar 
(or axion) piece of the magnetoelectric susceptibility of Cr203 violates parity and time inversion 
invariance and doesn't contribute to the electromagnetic energy. In Sec lVIIl we mention other sub- 
stances that, besides Cr203, carry an axion piece, and in See l Villi we will discuss the implications 
of our result to other disciplines within physics and electrical engineering. 

In this article, we will base our considerations on the 4-dimensional tensor analytical formalism 
as described by Post ll93ll 116 ill , which can be understood by experimentalists and theoreticians 
alike. Accordingly, the time coordinate = t has the dimension of 'time' whereas the spatial 
coordinates x^, x^ , x^ are related to the dimension of 'length' . Only intermittently we will mention 
the formalism of exterior differential forms, as it is used in the book l|26|] . The equations in our 
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paper are quantity equations that are valid for an arbitrary system of units. However, if in SecslTVl 
and|V]we turn to experiments, then we use SI and sometimes Gaussian units, which are widespread 
in the literature. 



n. LOCAL AND LINEAR MAGNETOELECTRIC MEDIA, THE PSEUDOSCALAR a 



A. Excitation and field strength 



Post [16 ih represents the 4-dimensional electromagnetic field tensors according to 

/ D2 Di\ 



(5) 



(6) 



-D, Hs -H2 
-D2 -H^ Hi 
\-D3 H2 ~Hi / 

^ — El — E2 — E^ ^ 

El —B2 

E2 — -B3 Bi 

\Es B2 -Bi / 

The coordinate indices ^u, z/, ... run from to 3. We have to distinguish carefully between the upper 
and lower indices for reasons of general 4-dimensional invariance. We use Post's conventions 
throughout, unless indicated otherwise. 

The field 0'^'^, the electromagnetic excitation, represents a tensor density of weight +1. E^„ is 
a tensor. In order to transform (25^*^ to a pseudo tensor with lower indices, we introduce (Einstein's 
summation convention assumed) 



G 



2 ^^iukX^ 



kX 



-G 



VfJi 1 



(7) 



with the totally antisymmetric Levi-Civita symbol l^^^x = 0, ±1. Here 'e^^^x is a pseudo ten- 
sor density of weight —1. We use e^^uKX for the dualization in ([7]), since this relation is invariant 
under proper as well as under improper transformations (i.e., those that include coordinate reflec- 
tions). We will denote pseudo tensors by a tilde throughout our article. By simple algebra we can 
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construct the corresponding matrix for G^,y: 

f Hi H2 \ 

-Hi Ds -D2 



G 



111/ 



-G 



i/fi 



(8) 



~H2 -D-i Di 
\-H3 D2 -Di / 

Whereas Ffj_^, the electromagnetic ^zeW strength, is a tensor, G^i, is a pseudo tensor, i.e., they 
behave differently under a reflection x° — x° or —x^ etc. In the language of dif- 

ferential forms, the electromagnetic field can be represented by two two-forms, the excitation 
G = G^^ dx^" A dx''/2 and the field strength F = F^^ dx^" A dx''/2, respectively. The G is a 
pseudo 2-form, also called twisted (or odd) 2-form, the F is an untwisted (or even) 2-form. 



B. The Maxwell equations 



The Maxwell equations in premetric form are [|61|] : 

d^F.x + d^Fx^ + dxF^, = . 



(9) 



If we introduce the dual of the electric current density, J^^x '■= ^^iuXk then the inhomogeneous 
equation can be transformed into the equivalent equation 



d^Gyx + duG\^ + d\G^y — J^ux ■ 
In a more condensed way, we may also write the Maxwell equations as 



d[iiGiyX] — -^J^iux , d[fj_Fj,x] — . 



(10) 



(11) 



where [fii^X] = — z^/iA + Xfiv =f • ■ ■ )' i-^-' ^ V^^^ (minus) sign occurs before even (odd) 

permutations. In the differential form language, the Maxwell equations (fTTI) read 



dG = J, dF = 0, 
with the 3-form of the electric current J = hJ^^x dx^ A dx" A dx^. 



(12) 



C. Constitutive relation 



The system of the 4 + 4 Maxwell equations Q for the 6 + 6 independent components of 
the electromagnetic field and Fxp is evidently underdetermined. To complete this system, a 
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constitutive relation of the form 

6^^^ = (B^'^iFxp) (13) 

has to be assumed. The constitutive relation (fT3] ) is independent of the Maxwell equations and 
its form can be determined by using experimental results. For vacuum, the constitutive relation 
("spacetime relation") is simple, 

= Yoy^F^" , (14) 

with g := det ^ and F^*^ := g^'^g"^ F^p. Here g^^^, are the co variant components of the 
metric of spacetime with signature ( — 1-+ +), its contravariant components g^^ can be determined 
via gp_x g^^ = 5^. For the excitation pseudo tensor, we have G^j^ = Y^'e^ycrr^f—g and, in 

exterior calculus, with the Hodge * operator, G = Yq*F. 

Turning to general magnetoelectric media, we assume with Tamm [|820, see also Post [16 ill , the 



most general local [l94ll and linear homogeneous relation 



0^- = ix'"'^-F.K, (15) 

where ■^^^"'^''^ is a constitutive tensor density of rank 4 and weight +1, with the dimension [x] = 
1/ resistance. Since both (S^^ and F^-k are antisymmetric in their indices, we have ■y^^'^"'^ = 
_^\uK(T _ _^uXcTK_ js^^ antisymmetric pair of indices corresponds, in four dimensions, to six 
independent components. Thus, the constitutive tensor can be considered as a 6 x 6 matrix with 
36 independent components. 

A 6 X 6 matrix can be decomposed in its tracefree symmetric part (20 independent components), 
its antisymmetric part (15 components), and its trace (1 component). On the level of this 
decomposition is reflected in 

36 = 20 © 15 © 1. 

The third part, the axion part, is totally antisymmetric and as such proportional to the Levi-Civita 
symbol, i^)-^^'^'^'^ ;= -^[Aj^ctk] _ ~~XuaK_ jj^g second part, the skewon part, is defined according to 
{2)^fj.uXp ._ i(^^tiuXp_^XptJ,u^_ If the constitutive equation can be derived from a Lagrangian, which 
is the case as long as only reversible processes are considered, then ('^)^^^'^'^'^ = o. The principal 
part fulfiUs the symmetries (i);^^'"^'^ = and = 0. The constitutive 

relation now reads 

= 1 ((1);^^^-'^ + (2)^A^-« + ae^'^''^) F,^ . (17) 

1 



In order to compare this with experiments, we have to split (fTTl) into time and space parts. As 
shown in [|27ll in detail, we can parametrize the principal part by the 6 permittivities e""^ = e^"". 



the 6 permeabilities fiab = f^ba, and the 8 magnetoelectric pieces 7'*;, (its trace vanishes, Yc = 0) 
and the skewon part by the 3 permittivities Ua, the 3 permeabilities m"-, and the 9 magnetoelectric 
pieces s^^. Then, the constitutive relation (fTTl ) can be rewritten as 

I^" = {e"' - e"'^ n,) E, + ( 7% + Sb" - E'' + 5 5" , (18) 

Ha = {fiat - eabcm') B' + {-^\ + sj ' 5^/) E, - a Ea . (19) 

Here e"*'^ = iabc = ±1, are the 3-dimensional Levi-Civita symbols. As can be seen from our 
derivation, 5 is a 4-dimensional pseudo (or axial) scalar, whereas Sc'^ is only a 3-dimensional scalar. 
The cross-term 7";, is related to the Fresnel-Fizeau effects. The skewon contributions m^, ric are 
responsible for electric and magnetic Faraday effects, respectively, whereas the skewon terms Sa 
describe optical activity. Equivalent constitutive relations were formulated by Serdyukov et al. 



117 511 . p. 86, and studied in quite some detail. 

According to Post [61], the pseudoscalar a should vanish for the vacuum and for all media. We 
will show in the next section that, in general, this is not the case. For Cr203 the pseudoscalar a 
turns out to be finite. 



III. THE ANTIFERROMAGNET CR2O3 AND THE THEORY OF 
DZYALOSHINSKH 



On the basis of neutron scattering data DSD and susceptibility measurements 115 ih of the an- 
tiferromagnetic chromium sesquioxide [l95ll Cr203, Dzyaloshinskii [|l5n was able to establish the 
magnetic symmetry class 3'm' of the Cr203 crystals. In accordance with these results, Dzyaloshin- 
skii developed, by starting from a thermodynamic potential quadratic and bilinear in E and H, as 
foreseen by Landau & Lifshitz [42], a theory for the electromagnetic constitutive relations for 
Cr203. We write them here as quantity equations that are valid in an arbitrary system of units: [|96|] 



Dnc = e±eoErc + axy/e^^H^ , 
Dy = e±eoEy + a_L^/e^^Hy , 
Dz = £\\ eoEz + a\\ y/EofioHz , 



(20) 
(21) 
(22) 
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and 



Bx = IJ'±IJ'oH^ + a±y/eoiJ,oEx , (23) 



By = fi±fioHy + a±y/eofioEy , (24) 
Bz = fJ^w fJ'oHz + a|| y/EofioEz . (25) 

The 2-axis is parallel to the optical axis of Cr203. Remember that y/eofio = 1/c and \/eo/Jto = 
1/Zq, with c as velocity of light and Zq as vacuum impedance. Here we have permittivities parallel 
and perpendicular to the z-axis of the crystal, namely e\\,£±, analogous permeabilities /U± and 
magnetoelectric moduli a\\,a±. Note that all these moduli are dimensionless, and this is true for 
all systems of units. Our dimensionless a's are different from the ones used by experimentalists 
and theoreticians up to now. We will discuss the transition to different systems of units in Sec lIV Al 
below. 

As we can see from ([5]) and ( [T?] ). we have to get (D, H) on the left hand side and (E, B) on 
the right hand side in order to end up with a constitutive law that is written in a relativistically 
covariant form. For this purpose, we resolve the last three equations with respect to H: 

Hx = Bx —J— Ex, (26) 

/i_L/iO /i± V /^O 

Hy = By —.j-^Ey, (27) 

/i_L/iO /i± V /^O 

Hz = Bz-^J—Ez. (28) 

/i||/iO /i| 

On substitution into dlQl) to we find. 




Dx = {e±-^] BoEx + —.[^Bx , (29) 
Dy = [e^-^] SoEy + —.[^By , (30) 




Dz = \ e\i-^] BoEz + ^J—Bz. (31) 

Now we have to compare with the local and linear constitutive relation (fTSl ) and (fT9l) . Since 
Dzyaloshinskii assumed that his constitutive relations can be derived from a Hamiltonian, it is 



clear that the skewon piece with its 15 independent components vanishes identically, see [|26|], 
Eq.(D.1.44). But this can be also read off from comparing ([H]) and with to dlB- The 
skewon pieces Uc and must be zero, since the Da in (|29l) to (|3TI) are proportional to Ea, and the 
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Ha in (|26|) to (|28]) are proportional to A similar consideration shows that sj' = 0, since Sa^ 
only provides off-diagonal pieces. Consequently, eqs. (fT8l) and (fT9l) reduce to 

(32) 
(33) 



Ha = ^-^B'-^\Eb-aEa, 



with 21 independent moduli. The permittivity matrix e""^ and the impermeability matrix are 
both symmetric and possess 6 independent components each, the magnetoelectric cross-term with 
which is tracefree, = 0, has 8 independent components. The 4-dimensional pseudo scalar 
(we call it also the axion parameter) represents 1 component. 

By comparing (|32l) and (|33l ) with the above equations (|29l) to (|3T|) and (|26l) to (|28l) , we can read 
off the permittivity 



£ — £0 
















^ 

£11 - 



and the impermeability 



(fll^ ^ 
f^l^ 

i^O /.-j 



For the magnetoelectric cross-terms, we have 

'y\B' + aB' 

and 

- j""^ E^-aE^ 
-7% Ey-aEy 
-7^ E,-aE, 




V /^o 

V /io ^ 



(34) 



(35) 



(36) 
(37) 
(38) 



(39) 
(40) 
(41) 
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Note that in the Cartesian coordinates used by Dzyaloshinskii we have 5^ = Bx etc., since the 
spatial metric is Euclidean with signature (+ + +). Thus, we are left with 

V /^o 

Y /io 

One of the triplets of equations (|36l) to (|38l) and (|39l ) to (|4T]) is redundant because of the vanishing 
of the skewon piece of 

The matrix 7 is traceless: 7^^. + 7^^ + Yz = 0. If we add up all three equations, we find for 
the pseudoscalar (or axion) piece 



7% + « 



(42) 
(43) 
(44) 



a 



a I a\\ 

/i_L fi\\ 




(45) 



Resubstituted into (l42l) . (|43]) . and (l44l) . the magnetoelectric 7 matrix becomes 



'^l 0^ 



7 6 



3 \fi± 




(46) 



7 



1 

that is, it has only nonvanishing diagonal components! The magnetoelectric matrix 7^*;, as well as 
the pseudo scalar a carry the dimension of 1/resistance. Conventionally, in the "magnetoelectric 
literature" the 7-matrix and a are collected in the "relativistic" matrix 



V f^o 



o^ 







MX 







0^ 

Ml 



(47) 



It is called relativistic, since it occurs in the context of the relativistic (E, B) system, see (|32|) . 

Since there are doubts in the literature about the correctness of Dzyaloshinskii's theory, see 
Lakhtakia |40], it is important to note that O'Dell S, PP- 115/116, and Janner [32], p.205, (see 



also Rado & Folen [165D and O'Dell 



) analyzed the crystal structure of Cr203 and determined 



the form the matrices e"-^, and '^^^a^b ought to possess. They found 



^ah , — 1 rel „ a 



(48) 
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where nonvanishing entries are denoted by • and *, respectively. A comparison with (l34l) . (l35l) . 
and (|46l) confirms Dzyaloshinskii's theory. 

Summing up: the nonvanishing magnetoelectric moduli for Cr203 can be determined with 
the help of eqs. (|45l) and (|46l ). Let us stress that £^^,e_i_, yU^, a^\, and a_i_, according to their 
definitions (|20l) to (|25l) . are measured in an external E and/or an external H field. 



IV. MAGNETOELECTRIC EXPERIMENTS WITH CR2O3 



A. Dimensions and units 



Since in the literature, which we need for extracting data, different systems of units are used, 
we want to underline again that usually all our equations are quantity equations, which are valid in 
all systems of units; only in the "experimental" SecsjIV]and|Vl we will turn to specific systems of 



units and some equations may be unit-dependent, see also Rivera [|69|1 in this context. We will go 
into some detail here, since these questions often lead to misunderstandings between theoreticians 
and experimentalists. 

A physical quantity Q is given by 

Q = {Q} [Q] ■ (49) 

Here {Q} is a numerical value and [Q] the physical dimension of the quantity Q. For instance, we 
have for a (1 -dimensional) displacement s, 

s = {s} [s] = {s} length, (50) 

where length is the dimension of s. So far, all of the equations in SecsjI]to|in]are quantity equa- 
tions. They interrelate "physical quantities" that consist of numbers and dimensions. Like in (|49l) 
and (l50l) . they are totally independent of any system of units. In many papers and books, the equa- 
tions are only valid in one system of units, they are numerical equations, like in Jackson's book on 



electrodynamics [31^, for example. Il97l1 In our paper, as in Post [I6III or in Il26ll . the equations are 
quantity equations and are always valid, independent of the units chosen. 

In a second step, if we relate our equations to measured values, we need a system of units. Then, 
for the example above, 

s = 15 m = 1500 cm «i 45 /eet = ... 

= {s}' m = {s}'" cm = {s}" feet = ... (51) 

12 



The numerical value {s} depends on the unit chosen. In fact, 

{s}' [s]" {s}' [s]" 



(52) 



{s}" H" {sV 

etc., that is, we have reciprocal proportionality. This is one of the fundamental laws of dimensional 
theory. The physical quantity s is invariant, i.e., it doesn't change, but its numerical value {s} may 
change according to the choice of the unit. 

In the center of our interest is the pseudoscalar 5. According to (|33] ), it has the dimension 

[Ha] current length 1 

\a\ = j—i = — X — — = — : . (53) 

[Ea\ length voltage resistance 

Moreover, the permittivity matrix e""^, according to (|32|) . has the dimension 

, ^ ^ charge ^ kngth ^ ^ I ^^^^ 

[Ea] area voltage velocity resistance 

-1 



and the impermeability matrix /i^^ , according to (|33l) . 

r -n [Ha] current area . 1 

A^nh = 77TT = -; 7- X — ; ■■ — = velocity x — . (55) 

[B"-] length voltage x time resistance 



In particular, we have 



[^o] = A- X , [l^o ^] = vel. X ^ 



vel. resist. resist. 
Accordingly, we can summarize these considerations in 




resist. 



(56) 



£"^1 \U.:\ 17"/, I \a\ 



1. (57) 



N [f^o^] W^o/l^o] IV^o/l^o] 
These quotients are dimensionless in all systems of units. Here Eq and jiQ are universal constants. 
They have in different systems of units different numerical values. We have y/eojio = 1/c and, in 
particular, 

^ 1 1 



/io Zq 376.73 n 



in SI, (58) 



' £2. = _^ in Gaussian units, 

Ho A-K 

' — = c in Heaviside-Lorentz ("rationalized Gaussian") units. 



Therefore, if you see and fiQ in our equations, like in (|2()1) to (1251) . it does not mean that we are 
in SI, but rather that we use quantity equations with correct dimensions. 
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The thermodynamic potential (an energy density) relevant for the magnetoelectrical effect con- 
tains cross terms between electric and magnetic fields. As a quantity equation, it reads, restricting 
ourselves to the linear regime. 



^?(E,H;T) 



t/e 



dim.— less 



(59) 



This linear part of a power series development is valid for "small" values of the Ea and Hb fields, at 
least relatively to the internal crystal fields. The a± and a\\ of Sec|lll]are components of the matrix 
a"''. Hence, strictly speaking, we should have put a star to all of them. However, for convenience 
we dropped these stars. 

We will concentrate here on SI and on the Gaussian system. In the Gaussian system of units, a 
mixed system consisting of electrical (electrostatic) and magnetic cgs-units, see, e.g., Sommerfeld 
[|8l|l . Panofsky & Phillips [|58|], or Jackson [|3l|], we have the following field redefinitions: 



^E,, = E,, , = 4nD„ 



-Hr. 



G 



•Ba = cBa 



(60) 



The speed of light c is instrumental for making the dimensions of '^Ea and '^Ba equal to each other, 
^Ea\ = ['^Ba] , and, similarly, ['^Da] = ['^Ha] . The An removes this factor from the Coulomb law. 



With the field redefinitions (1601) and with the convention in the Gaussian system 



ca 



ab 



Eq.(l59l) can be rewritten as (see Landau-Lifshitz [|42l] and Dzyaloshinskii [|l5|] for the Gaussian 
system) y 



-g{E,H-T) = ■■■ + V^'E,^'Hb = --- + ^o ''E^^'H, 



ST afe Slzr Sir 



t/e 

I 1 G^ab Gp Gtt 

An 

dim.— less 



dim.— less 



(61) 



In SI, fh"^] = s/m, f^Ea] = V/m, and [^^Hh] = A/m. Then, 



[SI, 



'9] 



VAs J kg 



rrfi rrfi ms'^ 



10 



9 



cm 



10 fg] . 



(62) 



On the other hand, in the Gaussian system, we have for the electric field 

l^^E^ 



cm m 



(63) 



and for the magnetic excitation 



An m An 
14 



(64) 




FIG. 1: The MEe effect (linear magnetoelectric effect with electric field-induced magnetization) of Cr203: 
Temperature dependence of the magnetoelectric components a\\ and a_i_ according to Astrov [3], see also 
the Tables 13] . 



Thus we have for the magnetoelectric moduli 

1 = c fa''''] 



3 X 10« - [V] 

s 



(65) 



Accordingly, we have the rule that multiplying ^^a"*, given in s/m, by c = 3 x 10® m/s yields the 
dimensionless Gaussian value '^a"*. Incidentally, in some papers Heaviside-Lorentz ("rationalized 
Gaussian") units are still in use, see, e.g., Borovik-Romanov & Grimmer [8], p. 139. 



B. Astrov, Rado & Folen, and Wiegelmann et al. 



In our task to determine the pseudoscalar a, we can take recourse to already published experi- 
mental data. Our main sources are Astrov Is] (Fig. l) for the electrically induced magnetoelectric 



effect (MEh), see also O'Dell [Svl] and JJ, 



ijleifj 



effect (called MEe in future) and Rado & Fole n 16411 fo r the magnetically induced magnetoelectric 



of Cr203 were used. In the MEe experiments [|3|, 



11 



2d 



63h . In both investigations single crystals 



6411 . Eqs.(|23l) to (|25l ) were verified {H switched 



off) and in the MEh experiments [|64ll Eqs.(|20l) to (l22l) {E switched off). In particular, Rado & 
Folen made both type of experiments and found that the magnetoelectric moduli a±_ and an for 
MEe experiments coincide with those of the MEh experiments. This proves the vanishing of the 
skewon part of the constitutive tensor -yA^'^'^ for Cr203. 



15 




75 100 125 150 175 200 225 250 275 300 325 

TEMPERATURE T, K 



FIG. 2: The MEh effect (linear magnetoelectric effect with magnetic field-induced polarization) of Cr203: 
Plot of a||(r) and a±{T) after digitalization and interpolation (small full triangles) of Fig.2.8, p.41, of 
Wiegelmann 18811 . (see also Fig.2, p. 143, Wiegelmann et al. [90]). The sign of a_L(T) was set negative 
according to Astrov t^. The curves (B -Splines) are only guides for the eyes. 



Accordingly, these experiments confirmed Dzyaloshinskii's theory for Cr203 below the spin- 
flop phase. Further experiments were then done mainly for the MEh case since (i) it is easier to 
conduct an experiment with very high magnetic fields rather than with high electric fields, (ii) 
Even at low i^qH fields, say below 1 tesla, the MEh effect needs no calibration of the measuring 
system, contrary to the MEe case! With the quasi-static method, as Da = —dg/dEa = a'^''Hb and 
Da = Q/Sa, a"^ = Q/ (SaHh). The charges Q, usually in the pC range, are measured with a high 
input impedance electrometer, the magnetic field /io-f^6 with a Hall probe and the area Sa of one of 
the electrode by taking a picture of it. (iii) Quasi-static as well as AC measurements can be done 
easily. For MEe quasi-static experiments, a SQUID must be used, see Kita |38, 3^. 

With either methods, MEh or MEe, (i) best results are obtained with gold electrodes rather than 
with silver paste ones, (ii) The resistivity of the sample must be high enough, especially when the 
MEh quasi-static method of measurements is used, (iii) One has to make sure that the antiferro- 
magnetic magnetoelectric crystal forms a single domain by cooling it in appropriate simultaneous 
electric and magnetic fields through the Neel temperature; this is the so-called magnetoelectric 



16 



annealing, see ll49ll . [|50h . and nSm . p. 124. Magnetoelectric crystals may also be (weakly) ferro- 
magnetic or ferrimagnetic. 

Our third main source of information are the measurements of Wiegelmann et al. [|90(], see also 



m 



8911 . He took magnetic fields B as high as 20 tesla and measured from liquid Helium up to 



room temperature. Wiegelmann et al. took a quasi-static magnetic field and thereby disproved 
explicitly claims by Lakhtakia [|40n that measurements with magnetic fields of some kilo hertz 
cannot be extrapolated to static measurements. The 0:22 values of Wiegelmann et al. [90] were in 
very good agreement with independent a^z measurement presented below (Sec lIVCl) . However, 
the sign of a±(T) relative to a\\(T) was left open. Hence we took that from Astrov [|30. 

The values of a^{T) and a\\ {T) of Fig. 2 are thus taken from Wiegelmann [|88D, Fig. 2. 8, p.41, 
(see also Wiegelmann et al. [90] (Fig. 2, p. 143)) after digitalization, interpolation, and correction 
for the relative sign. These values are given here in SI units in ps/m (picosecond/meter). 

Independently, see SeclTVl Rivera measured quasistatically at 133 Hz a\i of Cr203 between 
1.6 K and 305 K. He normalized that at T = 275 K, the temperature of the maximum value of 
a\\. He found the maximum value of a||(at 275 K) = 4.13 ps/m. As we saw in (|65l) . ^^a"'', given 
in s/m must be multiplied by the speed of light||99|] in order to yield the dimensionless ^^af'. 
Dropping again the star, we have 



a\\ (at 275 K) = 4.13 x 10"^^ — x 2.99792 x 10^ — ^ 1.238 x 10" 

m s 



(66) 



C. Unpublished measurements of Rivera (1993) 



One of us conducted quasi-static MEh experiments on Cr203 some years ago which were only 
briefly mentioned in [69]. 

After grinding and polishing a crystal, polarized light microscopy was used, see, e.g., 
Wahlstrom [|87n or Hartshorne [|24ll . to check the orientation of a perpendicular cut to the opti- 
cal axis. Actually, conoscopy allowed such tests (Fig. 3) in the near infrared, 260 /im thick Cr203 
crystals being absorbing in visible light. The symmetry was found to be uniaxial to the accuracy 
of the conoscopic method, consistent with the symmetry 3'm' [9], although now it seems to be 3' 
[|l4ll . The optical axis was computed to be inclined less than 3° away from the normal to the cut. 
Incidentally, this platelet was then sent to Wiegelmann in Grenoble, see [|90n . 

For the MEh experiments of ol\\{T) presented below, semi-transparent gold electrodes were 
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FIG. 3: Conoscopic pictures, in the near infrared, of a Cr203 platelet at room temperature and between 
linear crossed polarizers, prepared for measuring the a|| (T) coefficient. We can see one or two isochromatic 
curve(s), the circle(s), the cut being slightly inclined, less than 3° relatively to the optical axis, the z axis. 
At the center of the cross, the melatope emerges as the optical axis with a possible rotation of the E vector 
of the hght along this axis, evidenced by a lighter center (Zeiss objective 125 x /1. 30 oil P, and condenser 
with head and oil, numerical aperture = 1.30). 

evaporated on both sides of a platelet with area (one side) = 4.70 mm^ (thickness th^ = 
260 fim). On Fig.4, this Cr203 platelet is shown mounted with its thin gold wires (0 40 /im) on 
Cu wires and then on two coaxial 50 low noise cables #1 and #2, on a stainless steel sample 
holder. On the top left (forefront) of the picture, we see the temperature sensor, a calibrated (1.5 K 
to 300 K) Carbon Glass Resistor (CGR-2000, "Lake Shore", about 3 mm). This sample holder 
was then inserted in a copper can with He exchange gas, in an Helium bath cryostat. 

As already mentioned, before the measurements, the crystal was always cooled using the so- 
called magnetoelectric annealing with DC magnetic and electric fields, see Martin, Anderson, and 



Schmid 



50 



2M in this general context. It was easy to apply an electric field because the 
electrical resistivity of the crystal was very high. We measured the charges Q with a low noise 
electrometer (Keithley, 642-LNFR\) and H with a Hall probe. The H was produced by an old 
12-inch Varian "V 4012-3B" electromagnet with Varian "Mark I" magnetic field regulator. We 
measured at more than 130 temperature values. Typically, the record of Q was as follows (quasi- 
static MEh method): for about half a minute at zero field, then increasing linearly H with time, 
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' — e: 



FIG. 4: a) A Cr203 platelet (Sz = ^^.l^mm? ,thz = 260 jim) for a\\{T) measurements, connected to two 
low noise coaxial cables #1 and #2. One can see, by reflection on the platelet, the gold wire and the Cu 
wire #2 to the right. Black lines were added on the right of this unique Polaroid picture for clarity, b) A 
schematic drawing of the set-up depicted under a). 

from to 10 kOe in about 2.5 minutes, maintaining H for half a minute at maximum, decreasing 
H linearly with time to and then maintaining if at for half a minute. This procedure was 



6911 ■ The final curve looks 



used to cancel any shift of the base line, if any, see also Rivera 1681 
like a volcano with a flat top and straight sides, proving the linear character of the MEh effect 
for ail of Cr203, see Fig. 5 for T = 291.7 J^. Note the very good signal over signal + noise ratio. 
Remember that lOe = 1000/(47r) A/m ^ 79.6 A/m and H = IQkOe yields B = ^qH = IT. 
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FIG. 5: Quasi-static MEh experiment for Cr203, charge vs. time. The field H varies linearly in time, from 
to IQkOe and from lOfcCfe to 0, proving the linearity of the magnetoelectric effect in Cr203, here for 
a||(at292 K). 

As Dz = Q/Sz = a\\Hz, we compute directly a\\ = Q/{SzHz). 

In Fig. 6, we show the curve of a\\{T)), in SI units, i.e., ps/m. The curve, around T = K, 
has an elbow which could come from a possible crystallographic and/or magnetic phase transition. 
Further experiments should be done below 10 K to clarify this point. Closely below Tn, from 
293 K to 304 K, we plotted (not shown) ln[a|| (T) /af^^^] vs. ln[(TN - T) /Tn]. It is a straight line, 
the slope gives the exponent ^ 0.34 (thus about 1/3). Astrov [3], from the cut for measuring 
a^{T) by a MEe experiment, found an exponent 1/2, as expected according to the Landau 
theory. 

In Table 1, Wiegelmann et al. ^90^ compared the values of the temperature T at a;|| =0 (where 
the sign changes) and T at a™^"" that were obtained by different authors. This is better than giving 
an error on the results of a because with the MEh effect, the larger error, apart from the one 
on the area of the crystal (about 3% - 5%), comes from the uncertainty about the success in the 
magnetoelectric annealing, which could also be influenced by the quality of the crystal. 

In Fig.7, we display Q{9) for H = 5 kOe and 10 kOe, also at T = 292 K, showing very good 
parts of cosine curves. The angle 9 is measured between the optical axis (the z axis) and the H 
field direction. This demonstrates that a^a; and are very small or even null, at least, at that 
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FIG. 6: Magnetoelectric susceptibility a\\{T) for Cr203, obtained by quasi-static MEh experiments, H 
varying at each temperature from to 10 kOe and back to 0, as described in the text. The interpolation 
between the points was made by means of B-sphnes and is only a guide for the eyes. 




FIG. 7: MEh experiments for Cr203 at T = 292 K and constant H fields. By rotating the electromagnet, 
the measured charges (continuous lines) follow cosine curves. Superimposed are normahzed computed 
points (open circles at = 10 kOe, open stars at = 5 kOe). This proves that azx and azy are very small 
or null. At9 = 0°, the field H is parallel to the optical axis (z axis). 
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FIG. 8: Magnetoelectric susceptibility a||(T) for Cr203 obtained by low frequency (133 Hz) AC MEh 
experiments (hz ~ ISCferms)- Knowing T w 85.0 A' at a||(r) = from quasi-static experiments and the 
base line above Tn, the curve was corrected for a small drift. 

temperature. This again supports Dzyaloshinskii's theory, see (|47] ). After transformation from 
rectangular to polar coor dinate s, we could obtain a continuous curve cos^ 9, similar to the one of 



Fig. 3 given by Astrov Oil 111 0011 for the case a\\ at T = 103 K. 

Fig. 8 concerns the MEh effect at a low frequency of 133 Hz. It represents a dynamic measure 
of «||(T), with an AC field ~ ISCterms (root mean square) superimposed to a low DC field 
Hz ~ 580 Oe, just to maintain the domains close to Tn. The rate of heating from 4.2 K to 330 K 
was +2 K/min. As shown in Fig. 4, the crystal was mounted "floating" between the coaxial 
cables #1 and #2. As Q = CU, by measuring with a "lock-in" amplifier (SR 530) the voltage in- 
duced by the AC magnetic field hz we have a measure of the QiT) which is proportional to a\\(T). 
To a first approximation, the capacitance C of the crystal was supposed to be independent of tem- 
perature. If we superimposed the normalized curves at Tmax of a\\{T), measured quasi-statically 
and dynamically, from 70 K to 270 K, the superposition is very good. The Neel temperature found 
dynamically (308 K) is 3K above the one found quasi-statically (305 K). Probably the heating 
rate (+2 K/min) was too high. 

This concludes the presentation of our MEh experiments on Cr203 which are in very good 



agreement with the ones presented, in particular, in ll90ll . 
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V. EXTRACTING THE PSEUDOSCALAR 5 



A. The permeabilities ij.± and /i 



The magnetic susceptibility tensor of Cr203, namely x := 1 — (1//^) or , in c omponents, Xab = 
5ab — l^ab^ determined by Foner [21] as a function of the temperature . [| 1 ill From the caption 
of his Fig. 8, we can take the static susceptibility perpendicular to the c axis xl of Cr203. The 
superscript g stands for specific (or mass) magnetic susceptibility. Foner used the old unit emu/ g. 
From the inside of the front cover of Landolt-Bornstein [43] we learn that 1 emu/ g = 1 cm^/g. 
Accordingly, 



Xi(at4.2i\:) = 2.24 x 10" 



cm 



9 



(67) 



Let us now determine the volume susceptibilit y y^. Again from Landolt-Bornstein we take = 
4:71 The de nsity for the mineral eskolaite. il 10211 containing 94% Cr203, is 5.23 g/cm^ and for 
Cr203 ceramics [|l03ll a bit less, namely 5.21 g/cm^. Then, in SI, 



Xl(at4.2ir) = Attx 2.24 x 10"^ x 5.22 ^ 1.47 x 10^-^ . 



(68) 



We can read off from Fig. 8 of Foner [211] that both, the parallel and the perpendicular suscepti- 
bilities at the Neel temperature Tn are about 13% higher than x]_(at 4.2 i^'). Consequently we 
find 



Xl(at Tn) ^ xUat Tn) ^ 1.13 x xl(at 4.2 K) ^ 1.62 x 10 



-3 



(69) 



O'Dell |157|], App.l, found the slightly higher value of ^ 1.64 x 10"^. 

Now we can determine the permeabilities: Below and close to the Neel temperature, we have 

1 



1 - r (at 



1 + X''(at Tn) ^ 1.00162. 



(70) 



At 4.2 K, we find 



/ix(at4.2ir) ^ 1 +xl(at4.2ir) ^ 1.00147. 



(71) 



Since x\\ is ev en smaller, we have /i ~ 1. Because Cr203 is antiferromagnetic, this was to be 
expected. [11041] 
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FIG. 9: The pseudoscalar or axion piece a, see Eq.©, of the constitutive tensor ^^'^'^'^ of Cr203 in units of 
Yq = I/Zq as a. function of the temperature T in kelvin; here Zq is the vacuum impedance which, in SI, is 
~ 377 ohm. 

B. The pseudoscalar (or axion piece) 5 

Now we can come back to (|45] ). Since ~ ~ 1, the pseudoscalar (or axion piece) of the 
constitutive tensor becomes at^^'^'^, with 

5~ 1 (2ax + a||) ./^. (72) 



It is the arithmetic mean of the trace of the magnetoelectric tensor '''^'o;'^;, of (|47]) . Going back to 
Fig. 2, we can take the values of a± and a\ \ and compute 5. The result is plotted in Fig.9. As we can 
see, for temperatures of up to about 163 K, the pseudoscalar is negative, for higher temperatures 
positive until it vanishes at the Neel temperature of about 308 K. For example, at 285 K, we find 
5max ~ 1.035 ps/m. Multiplied by c we get 



5max (at 285 K) ^ 3.10 x 10"* J — 



4 / ^0 _ , 1 



/io 3226 







ST 11 

8.22 X 10-^ - ^ . (73) 

1.216 Ml] ^ ^ 
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VI. PSEUDOSCALAR OR AXION PIECE OF THE MAGNETOELECTRIC SUSCEPTIBILITY 
OF CR2O3 VIOLATES P AND T, IT CARRIES NO ELECTROMAGNETIC ENERGY DENSITY 



As has been pointed out by Janner [|33ll . amongst others, the violation of the invariances under 
space reflection (parity P) and time inversion (T) of the corresponding crystal, in our case Cr203, 
are necessary conditions for the emergence of the magnetoelectric effect. The same is true for the 
emergence of the pseudoscalar or axion piece a. The constitutive relation for the axion piece alone 
can be read off from (l32l) and (133]) as 



Ha = —a Ea , 



(74) 
(75) 



see 1261], Eqs.(D.l.l 12) and (D. 1.1 1 1). If we denote, as in crystallography, a space reflection by 1 
and a time inversion by 1', see Janner [|33n . then we have 



I'D'' = D'' , I'Ha = -Ha , I'Ea = Ea , 



-B" 



(76) 
(77) 



see also Marmo et al. (48']. If we now apply a space reflection to (l74l ) and (l75l) , then they transform 
into their negatives. 



D"" = -aB", 
Ha = +aEa. 



(78) 
(79) 



The same is true for a time inversion. In other words, the constitutive relations for the axion piece 



violate P and T invariance, see also [|l6l . 



85h . They are only invariant under the combined PT 



transformation. Consequently, the violation of P and T invariance is an essential characteristics 
of the axion piece 5. 

The first, to our knowledge, who tried to utilize the constitutive laws (l74l) and (l75l) for vacuum 
electrodynamics was Schrodinger [|74|] . p.25. He made the ansatz, compare ([T71) . 

1 







yXuCTK 77' 



(80) 



A look at dH]) and ^ shows that this, for a = \, yields the laws (l74l) and (l75l) . But Schrodinger 
rejected it as being unphysical if taken for vacuum electrodynamics. 
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Another property is characteristic for the axion piece: It doesn't contribute to the electro- 
magnetic energy density. This can be seen easily since the energy density in electrodynamics 
is I (D^Ea + HaB"-). If (1741 ) and (1751) are substituted, this expression vanishes. But even more 
so, also the energy flux density vanishes. In order to prove this, we turn to the energy-momentum 
tensor Ix" of the electromagnetic field that is built up from the energy density To°, the energy flux 
density Tq^ (with 6 = 1,2, 3), the momentum density Ij^ (with a = 1,2, 3), and the momentum 
flux density according to 



1 



V 



energy d. energy flux d. 
momentum d. momentum flux d. 



(81) 



To^ is also called the Poynting flux and %a the Maxwell stress. The energy-momentum tensor 



reads [cf. Post [16l|], eq.(9.55)] 

"Ix' = Zbl-^''''Fx., (82) 



with 



£ := ^ ©'^^F.. = \{H^B- - D^E,) ; (83) 



the last equation can be read off directly from ([5]) and Q. With some algebra, (l82l) can be rewritten 

as 

Xa"^ = \e''^'''{Gx^Fp, - Fx^G,^) , (84) 



4 



see Pa] . eq.(B.5.40). Since the laws (1741) and (fTSl) can be put together in the manifestly covariant 
form 

G^, = aF^,, (85) 



we see immediately from (1841) that 

TA''(of axion piece 3) = . (86) 
Thus, in particular, the electromagnetic energy density To° of the axion piece vanishes. 



VII. OTHER SUBSTANCES AND SYMMETRIES PERMITTING MAGNETOELECTRICITY 
WITH THE AXION PIECE 

In the present article the relativistic analysis is based on data of the antiferromagnet Cr203, 
because it represents so far probably the best studied magnetoelectric material and has diagonal 
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components of the linear magnetoelectric effect tensor a, see (l47l) . However, other materials and 
symmetries could have served the same purpose, in principle. Among the 122 Heesch-Shubnikov 
point groups 58 ones are permitting the linear magnetoelectric effect Ivoj], and therefrom 32 ones 
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720. Strictly speaking, our 



possess diagonal components of the magnetoelectric tensor a [Si, 
magnetoelectric tensor '^'^'a in (1471 ) belongs to the EB scheme, see (|32l ) and (|33l) . whereas the 
corresponding tensor in the literature [|73n is the one of the EH scheme. However, as we can see 
in (l47l) . because of ^ 1 the differences are marginal and don't touch our arguments. 

One can di sting uish three types of diagonals (for the complete set of magnetoelectric tensors 
see, e.g., refs. |69], 72], and [8], Tablel.5.8.1, for the examples cited, see ref. Isl], Table 1.5.8.2, 
except for [|23ll and 



1. ) 19 point groups with: an ^ ^22 7^ 0:33 

Examples: 

Point group m'm'm' : DyA103, GdAlOs, TbA103 
Point group m' : ^367013! 
Point group m'm'2 : CU3B7O13CI 

2. ) 8 point groups with: an = 022 ^ 033 ("11 = "22 = a±, "33 = a||) 

Examples: 



Point group 3'm': [CrsOg,] [11051] Nb2Mn409, Nb2Co409, Ta2Mn409, TaaCoiOg 



Point group 3': Cr203 [|230 



Point group A/ m'm'm': Fe2Te06 

3.) 5 point groups with: an = 0:22 = ^33 
Examples: 



Point group 4'3m' (expected): Cr3B70i3Br, Cr3B70i3l UM 



Thus, it is clear that the pseudoscalar a occurs in quite a number of different substances. Its 
existence can no longer be denied. 
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VIII. DISCUSSION 



The structure of the constitutive law ([85]) is not unprecedented, as we already discussed in [|27h . 
In electrical engineering, in the theory linear networks, more specifically in the theory of two ports 



(or four poles), Tellegen |[83, 



8411 came up with the new structure of a gyrator, which is defined via 



Vi = -SZ2 , 

V2= sii, 



(87) 



where v are voltages and i currents of the ports 1 and 2, respectively. Let us quote from Tellegen 
ll84ll . p. 189: "The ideal gyrator has the property of 'gyrating' a current into a voltage, and vice 
versa. The coefficient s, which has the dimension of a resistance, we call the gyration resistance; 
1/s we call the gyration conductance." The gyrator is a nonreciprocal network element. 

If we turn to the electromagnetic field, then because of dimensional reasons the quantities 
related to the currents ii, 12 are the excitations D'*, Ha and the quantities related to the voltages 
f 1, V2 the field strengths Ea, -B". Then we find without problems straightforwardly the relations 



Ea = —S Ha, 

3"= sD". 



(88) 



If we rename the admittance s according to s = 1/5, then (f87l) and (l74l),(|75]) coincide. Without the 
least doubt, the gyrator is in the theory of two ports what the axion piece is in magnetoelectricity. 
The axion piece 'rotates' the excitations, modulo an admittance, into the field strengths, as the 
gyrator the currents into voltages. 



These analogies or rather isomorphisms carry even further. In 2005, Lindell & Sihvola [|45 . 
46D, see also [|44fl. introduced the new concept of a perfect electromagnetic conductor (PEMC). 
It also obeys the constitutive law G^j^ = a F^,^ or (|88l) . The PEMC is a generalization of the 
perfect electric and the perfect magnetic conductor. In this sense, it is the 'ideal' electromagnetic 
conductor that can be hopefully built by means of a suitable metamaterial, see Sihvola 1770 ■ The 
pseudoscalar a is called Tellegen parameter by Lindell et al., see ll47ll . p. 13 (for a more general 
view, see BTSD); artificial Tellegen material has been produced and positively tested by Tretyakov 
et al. 1860 . amongst others. 



Continuing with our search for isomorphisms, we turn to axion electrodynamics, see Ni [|53ll . 

ItinHH. 



Wilczek [.91.1 . and, for more recent work. 



If for vacuum electrodynamics we add to 
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the usual Maxwell-Lorentz expression specified in (fT4l) an axion piec e pat terned after the last term 
in (fTTI) . then we have the constitutive law for axion electrodynamics. [| 1 0611 



= — 



(89) 



Alternatively, with the excitation pseudotensor (|7]) we find 

1 



G 



2Zr 



and, in exterior calculus. 



G= 1^ +5 



(90) 



(91) 



We discussed this 'spacetime relation' and also the corresponding Lagrangian in some detail in 
[|26|l and [|27|]. The Hodge star operator * is odd; it transforms a form into a twisted form, and 
vice versa. Therefore we could also denote it by * (we don't!). In Cr203 we had a ^ 10^'^ /Zq. 
It is everybody's guess what it could be for the physical vacuum. In elementary particle theory 
one adds in the corresponding Lagrangian also kinetic terms of the axion a la ~ g^^^d^a d^a and 
possibly a massive term ~ 5^. However, this hypothetical P odd and T odd particle has not 
been found so far, in spite of considerable experimental efforts, see Davis et al. 11211 and references 
given. 

The axion shares its P odd and T odd properties with the a piece of Cr203, with the gyrator, 
and with the PEMC. One may speculate whether an axion detector made of Cr203 crystals could 
enhance the probability of finding axions. 
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